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Abstract 



We use an instanton approximation to the continuous-time spin coherent- 
state path integral to obtain the tunnel splitting of classically degenerate 
ground states. We show that provided the fluctuation determinant is care- 
fully evaluated, the path integral expression is accurate to order 0{l/j). We 
apply the method to the LMG model and to the molecular magnet Feg in a 
transverse field. 



This paper is dedicated to the memory of Victor Belinicher, who was lost 
when Siberia Airlines flight 1812 was shot down over the Black Sea, Oct 4th 
2001. Victor made many contributions to physics, in particular to the spin 
tunneling problem. 
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I. INTRODUCTION 



One of the most convincing demonstrations of quantum effects in a near-macroscopic 
system is provided by the accurate measurement [|I| of the level splittings in the molecular 
magnet Fes- These splittings are caused by the quantum tunnelling of the direction of the 
large (J = 10) molecular spin between two classically degenerate energy minima. 

The natural tool for studying such spin tunnelling should be the spin {SU{2)) coherent- 
state path integral. It is easy to establish that this formalism gives a good qualitative 
description of the tunnelling process — including the dramatic topological quenching of 
the tunnelling ^ that makes the Fes results so interesting. Unfortunately, a straightforward 
application of the spin coherent-state path integral to compute the semiclassical propagator 
or the tunnel splitting yields results that are incorrect beyond the leading exponential 
order. 

Although there do exist other path integral approaches which find the splitting correctly 
PJ^, the resulting calculations tend to be intricate, and the simplicity seen in the conven- 
tional Schrodinger particle case is lost. These problems have lead to the spin coherent state 
path integral acquiring a reputation for being mathematically ill defined — or at least harder 
to deal with than the conventional Feynman path integral, whose mathematical subtleties 
have been well studied. 

Recently, however, it has begun to be appreciated that the problem with the spin coherent 
state calculation is simply that the fluctuation determinant has an "anomaly" , and that, once 
the "extra phase" provided by the anomaly is taken into account, the coherent state path 
integral gives correct answers. This extra phase seems to have been originally discovered in 
the 1980's by Solari UTO), but the significance of his result was not widely appreciated. It 
was then rediscovered by one of the present authors [|lT| and also by Vieira and Sacramento 



12|. The interpretation of the extra phase as an "anomaly" is due to the remaining authors 



of the present paper [T^ . 

These previous discussions of the extra phase were restricted to the case of quantum 
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evolution between generic values of the classical degrees of freedom. However, when we 
calculate tunnel splitting, the endpoints of the instanton path lie at local minima of the 
classical energy and, just as in the Schrodinger particle case, the Jacobi fluctuation operator 
has a zero mode which makes the inverse of its determinant singular and the general formula 
for the propagator inapplicable. Thus our earlier work had not fully established that the 
spin coherent state path integral gave the correct result for the tunnel splitting. This we do 
in the present paper. 

In the next section we provide a brief review of the spin coherent-state path integral, 
including the correction to the fluctuation determinant prefactor. In sections three and four 
we discuss the complications that ensue when there is a zero mode and provide a general 
formula for the one-instanton contribution to the tunnelling amplitude. In section five we 
apply this formula to the relatively simple case of the Lipkin-Meshkov-Glick (LMG) model 
and in section six we evaluate the tunnel splitting for a realistic model of Feg. 



II. SPIN COHERENT STATES 



We follow the conventions in [|l^ and define our spin coherent states to be 



\z)=exp{zJ+)\j,-j), (2.1) 

where |j, —j) is the lowest spin state in the 2j + 1 dimensional representation of SU{2) and 
J+ is the spin algebra ladder operator obeying 

J+\j,m) = ^Jj{j + l)-m{m + l)\j,m + l). (2.2) 

The variable 2; is a stereographic coordinate on the unit sphere with 2; = at the south pole 
(spin down direction) and z = 00 at the north pole (spin up). 

These coherent states are not normalized, but depend holomorphicly on z. This means 
that matrix elements such as {z'\0\z) are holomorphic functions of the variable z, and anti- 
holomorphic functions of the variable z'. 
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The inner product of two coherent states is 

{z'\z) = {l + z'zf^, (2.3) 



and they satisfy the overcompleteness relation 

2j + 1 f (Pz 



\z){z\. (2.4) 



TT J (l + Zz)2i+2 

Here d?z is shorthand for dxdy. The factor 1/(1 + ^z)^ combines with this to make the 
invariant measure on the the two-sphere. The remaining factor in the integration measure, 
1/(1 +^z)2-', serves to normahze the coherent states. 

We may use the overcompleteness relation to derive a formal continuous-time path inte- 
gral representation for the propagator 

K{ljXuT) = {Q\e~'^^\Q. (2.5) 



We insert intermediate overcompleteness relations into ( |2.5| ) and consider the limit N 
oo. This leads to the path integration formula ||rT| 



K{Cf. a T) = df^iz, z) exp{5(^(t), z{t))}, (2.6) 



where the path measure d^ is 



and the action S{z{t), z{t)) is 

S{-z{t), z{t)) = J {ln(l + Cfz{T)) + ln(l + z(O)O)} + - ^H{z, ^)| dt. (2.8) 

The c-number Hamiltonian, H{'z,z), is obtained from the operator H by 

H{z,z) = {z\H\z)/{z\z). (2.9) 

The paths z{t), z(t) obey the boundary conditions z{0) = Q, z{T) = (j, but ^(0), z{T), 
being actually z{0 + e) and z{T — e), are unconstrained, and are to be integrated over fl^ . 



The manipulations leading to the continuous time path integral are heuristic, but with 
careful treatment the formal path integral should be as useful as the familiar configuration 
space Feynman path integral. In particular the semiclassical, or large j, propagator can be 
obtained from a stationary phase approximation to the path integral p3|] . 

The stationary phase approximation requires us to seek "classical" trajectories for which 
S remains stationary as we vary the functions z{t) and z{t). These stationary paths will 
generally be complex. If we write z as x + iy and 1 = x — iy, then, except in special cases, X 
and y are not real numbers. In particular there is no requirement that ^(0) be the complex 
conjugate of z(0) = Ci, nor that z{T) be the complex conjugate of 'z{T) = C^j. Bearing this in 
mind, we make variations about a chosen path, and keep track of all boundary contributions 
resulting from integrations by parts. We find that 

l + (.,z(T) 1 + 2(0)C, 

- r [jM^ - ) - (- - ) ) 

Demanding that this change in the action be zero requires the trajectory to obey Hamilton's 
equations 

...,:(l±2£)!fl, .^-><i±H)!^, (,n, 

2j oz 2j oz 

together with the conditions 6z{0) = and 5^(T) = 0. We can therefore impose the 
boundary conditions z{0) = Q, z(T) = (j^, but ^(0) and z(T) are free to vary, and so are 
determined by the equations of motion. This is important because Hamilton's equations are 
first order in time and we cannot simultaneously impose initial and final conditions on their 
solutions. 

The dynamically determined endpoints can also be read off from the Hamilton- Jacobi 
relations that follow from ( |2.1U ). These are 



dS^ ^ 2jz{T) dS^ ^ 2jz(0) 

The Hamilton- Jacobi relations also tell us that 
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OSn] dSc] , 

showing that is a holomorphic function of Ci, and an anti-holomorphic function of C/. 
These analyticity properties of S^x coincide with those of K. This is reasonable since exp S^y 
is the leading approximation to if, and we would expect analyticity to be preserved term- 
by-term in the large j expansion. Finally, we have the Hamilton- Jacobi equation 

^ = -zii(C;,;^(T)). (2.14) 



In 13 1 we showed that after we compute the Gaussian integral over small fluctuations 



about the stationary phase path the resulting semiclassical approximation to the propagator 



IS 



^scl(C/,0,T)= 1 'd0^) exp|5'el(C/,0,r) + -y^ 0SK(t)t^tj 



(2.15) 



or a sum of such terms over a set of contributing classical paths. In this expression 



is the "extra-phase" discovered by Solari, Kochetov, and Vieira and Sacramento. 

The form ( p.l5|) is valid only if the prefactor is finite. When we compute instanton 



contributions to tunnelling there is a zero mode in the quadratic form for small fluctuations, 
and the resulting divergent integral over this mode is to be replaced by an integral over a 
collective coordinate labeling the instant that the tunnelling event occurred. This we will 
describe in the next section. 

III. DEALING WITH THE ZERO MODE 

As is usual in calculating tunnelling effects, it is convenient to perform the computations 
in Euclidean (imaginary) time. For the sake of symmetry we will take the time evolution as 
running from — T/2 to T/2 and the propagator ( p.l5|) becomes 



1 f /• 

K{Cf,Q,T) = [D{T)]-2 exp IS,, + \ J 



T/2 

I 

T/2 



dr 



(3.1) 



where again 0sk is the integrand of the Solari-Kochetov phase 

lfd(l + zzfdH dil + zzVdH' 



^SK 



(3.2) 



2 \dz 2j dz dz 2j dz J 
evaluated along Zci{t), Zc\{t), and D{T) is the fluctuation determinant. The latter may be 



found by the "shooting method". As explained in [|13|, this involves solving the equation 



B{r) -dr + A{r) 
dr + A{t) B{t) 



(3.3) 



where 



A 



9SK 

B 
B 



lfd{l + zzfdH d{l + zzfdH" 
2\dz 2j dz dz 2j dz 



d (1 + 



zzY 



dH 



dz 2j dz ' 
d (1 + zzfdH 



(3.4) 



-T/2) 



'0' 



(3.5) 



dz 2j dz ' 

with the initial condition 

,'ipL/ -T/2 

Given the solution of this equation, we read off the determinant as D{T) = iIji^{T/2). In real 
time, and when there are no problems with zero-modes, this recipe leads to the prefactor 
appearing in ( |2.15| ). 

Now assume that the coherent states \zi) and \zf) represent spins pointing along the 
directions of two equal-energy global minima of the Hamiltonian H. Because the gradient 
of the energy vanishes at both ends, the classical path joining Zi to Zf has the character of 
an instanton: as the total time taken to traverse the path becomes longer and longer most 
of the motion still takes place in an "instant", a flxed period short in duration compared 
to the total. When T becomes inflnite, the epoch of this "instant" is arbitrary and this 
leads to a zero-eigenvalue mode in the fluctuation operator. Thus D{T) is formally zero. 
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The problem of dividing by the square root of zero is avoided by introducing a collective 
coordinate for the tunnelling epoch, and the formal infinity in the one-instanton contribution 
to the propagator becomes a factor of T. 

The classical instanton solution can be written Zci{t — tq), Zci{t — tq) where tq is the epoch 
at which the tunnelling occurs. Since, in the large T limit, the action for the tunnelling event 
is independent of tq, the normalized zero mode is 

^0 = _ = . , (3.6) 

where g is chosen to make 

/ mim,dT= Ul + ,ijl\dr = i. (3.7) 

J-TI2 J-T/2 ^ ' 



The divergent Gaussian integration over the coefficient of the zero mode is replaced by 
an integral over possible tunnelling epochs tq by inserting a factor of 

\ fT/2 [ dT\ 1 
1 = -= / rfro — exp -— ^'(ro) (3.8) 

into the path integral, with the choice 

Hro) = / dr'-— -%{t' - tq) , (3.9) 

J-T/2 1 + 2;c1^ci(t' - To) \z{t')J 

and then proceeding in a manner similar to that used for quantum mechanical instantons in 
the Feynman path integral p!6| , p!7| : we first set z = Zci{t — tq) + 6z{t — Tq) and similarly z. 
Next, after observing that everything depends only on the combination r — tq, we change 
variables r — tq — r. The integral over tq is then trivial and gives a factor of T. Meanwhile, 
after an integration by parts and ignoring the fuctuations of {z,z) about (-2ci,^ci) which are 
of higher order, the Jacobian factor becomes 

^ = / dr'^—^ = — . (3.10) 

OTo J-T/2 l + ^ci^^ci yt^i(r')/ V9 

The quadratic term in the exponent is a projector onto the zero mode and replaces the 
vanishing eigenvalue by l/2ja. The net result is the replacement 



[D{T)\ 



T 




D{T) 
An 



(3.11) 



where Ao(T) is the eigenvalue that vanishes as T becomes large. 

The desired ratio, Det'(I/) = D{T)/\q, is equal to ip j;^{T /2) / Xq{T). We will not have to 
obtain '?/'^(T/2) and Ao(T) separately. 



The eigenvalue problem is 



A^; 



(3.12) 



where L is the same operator as in ( |3.3| ), but with boundary conditions ipx{—T/2) = 
?a(T/2) = 0. 

For finite T the shooting method solution, "^l, is close to, but not quite equal to, the 
"small-eigenvalue" eigenfunction, \I'ao- Although obeys the boundary condition at r = 
— T/2, it does not quite obey the boundary condition at r = +T/2. In turn "^Xo is close to, 
but not quite equal to, the infinite— T zero-eigenvalue mode, \E'o- 

Now \E'o obeys the equation I/\E'o = 0, but no particular boundary conditions at ±T/2. 
There is a second solution of this equation, Sq = {C,o,C,oY- The Wronskian of these solutions 

Mr) Ur) 
Mr) Ur) 

is independent of r. Next we observe that the differential equation ( 3.12| ) can be converted 
to an integral equation 



(3.13) 



^A(r) = ^L(r) + A / dr'Gir, r')^A(r'), 

J-T/2 

= ^l{t) + ^ f_^^^dr' [^o(r)H*(r') - So(r)M/* (r')] ^x{r'). 



(3.14) 



Since hir) obeys the boundary conditions at —T/2, and the integral vanishes at this point, 
we can find the eigenvalues A by requiring that the lower component at of ^\ vanishes at 
r = T/2. We are only interested in solutions where A = Aq is very small. Because of this we 
can approximate the '^x{t') appearing in the integral in ( p. 14] ) by the zeroth-order solution. 
In this way we see that 
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1 fT/2 



dr ^o(T/2)4(r) - ^iT I2)%{r)\ ^/^(r). (3.15) 

-T/2 L J 



Ao(T) y-T/2 

The integral in ( p.l5|) may be evaluated using only the asymptotic behaviour of \E'o and 



Ho, which involve z^i and Zci- This asymptotic behaviour depends only on the form of the 
Hamiltonian in the neighbourhood of the endpoints. 

In all cases we consider the instanton solutions have the property that ^ci = ^ci their 
endpoints. Here the asterisk denotes a true complex conjugate as opposed to the formal 
conjugate denoted by the bar. The coincidence of the formal and true conjugate occurs 
because these endpoints lie on the real unit sphereQ Taking this observation into account, 
we parametrize the Hamiltonian in the vicinity of the initial stationary point in terms of 
two frequencies, 00x^2-, as 

2i 



;i + ztz. 



\2 



ui{z ~ Zi){z - z*) + luj2{z - ZiY + luj;{z - z*f . (3.16) 



Since H{z, z) is real, so is uJi. Also, because the initial point is an energy minimum, we must 
have Ui > \uj2\- We can therefore define a real mean frequency, u, by 

A similar expression holds at Zf with the same values of uJi and 002 provided the degeneracy 
in the Hamiltonian is due to some symmetry. (There might be an extra phase factor in uj2, 
but this makes no difference to the subsequent calculation). 

As r becomes large and negative, B — > U2, B ^ uj^ and A = (psK ^^i, so we see that 

>o \ / ^0- \ / ^0 \ / ^0- 

e"^; _ H - (3.18) 



where 



^This is a not a trivial observation: in a more accurate representation of the Feg problem which 
includes fourth order anisotropy terms there are additional instantons for which this fails to be 
true. 
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UJ2 —UJ + UJl 



>0- 
>0- 



0. 



(3.19) 



There is an analogous relation for (^o-;'Co-)*- Wronskian to connect \E'o- 

with So-, so everything can be expressed in terms of W and the normalization g. Similar 
remarks apply to \l'o+ and Ho+. If we write 









< 







and apply the boundary condition at —T/2 we can find a and /?, and hence 

Inserting this into ( |3.15| ) and noting that the ipo ipQ terms dominate, we find 

1 



Ao(T) 



or, 



2 g^T 



Ao(T) 7/;o-V^o+4c^'' 
Thus the one-instanton contribution to the propagator is 



T/2 

K{zf, Zi, T) = exp \Sd+ 2 I 4>SKdT 




0+ 



\UJ2\ 



{2uTe~2 



iuT\ 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



Note that tpo, ipQ are proportional to y^, thus ^/g drops out and we can simply put g = 1 
in the sequel. 
Let 



then 



— ^ , ~P — UJ' 



-oo 



0+ 



T +00, 



(3.25) 



(3.26) 
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with 



N ^{l + ZiZi){l + ZfZf). 



(3.27) 



Using this we can write 



r r 

K{zf,Zi,T)^explS^i + ^ J 




C+C- 



1 



(3.28) 



IV. EXTRACTING THE ENERGY SPLITTING 

Again assume that the coherent states \zi) and \zf) represent spins pointing along the 
directions of two equal energy global minima of the Hamiltonian H. Let be the 

approximate (tunnelling-ignored) energy eigenstates localized near these minima. These 
should have their phases chosen so that when tunnelling is included the eigenstates become 
the linear combinations 

\i^±)^^m±m). (4.1) 

If the energies of these states are 

E± = E,, ± iA, (4.2) 
and Oq, = {za\ipa), then as T becomes large the coherent-state propagator 

K{zf,Zi,T)^{zf\e-^^\zi), (4.3) 

is given by 

K{zf,Zi,T) ^ a/a*e-^^-^sinh(|Ar), 

/i 1 A3t3 \ 

-aKe---"(^AT+i^ + ...j. (4.4) 

We will find the energy splitting, A, by evaluating K in the one-instanton approximation 
and comparing with this expression. 
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It is necessary to find expressions for the ampfitudes Oj and a/. These are obtained by 
looking at 

K, = {Zf\e-^^\z,) ^ |a,re-^-^, (4.5) 

and 

Ki = {z^\e-^^\zi) ^ la^l'e-^-^, (4.6) 

both evaluated in the harmonic approximation. This evaluation is performed in the ap- 
pendix. This results in 



i^/ = (l+^/-/)^^/-^e-i(™)- (4.7) 



and a similar expression for Ki. Thus 



e5ci+|/T52(<^sK-^i)rfr r— ^ 
Now 

2uj{ijj -\- uji) 2uj 



C+C- 



<^2 



2 



so finally 



where 



and 



(4.8) 



(4.9) 



A = 2uj^fPe\ (4.10) 



P = 4— C.C- (4.11) 



/OO /"OO 

a^,(T) + I / (0SK - a;i)rfr, (4.12) 
-OO J —OO 

where a^z is the kinetic term 

awz(T) = — (4.13) 

i + Zc\Zc\ 

in the classical action — the boundary terms having cancelled with the {l+'ZfZfYil+'ZiZiy 
in the denominator. 
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V. THE LMG MODEL. 



In this section we will evaluate the tunnel splitting in the relatively simple case of the 
Lipkin-Meshkov-Glick (LMG) model ll| . 



We will take the LMG Hamiltonian to be 

H = 701 + J') + ^, (5.1) 

with w > 0. Since + = 2(J^ — Jy), we see that the classical minima lie along ±y. The 
Hamiltonian which appears in the path integral is 

(z\H\z) r- z'^ + z'^ jw 
H{z, z) = \' , '/ = V2jw- + 5.2 

By setting dH/dz = dH/dz = 0, the classical minima are found to be at the points 

{z,z) = (5.3) 

which correspond to the ±y directions of the Cartesian axes. The explicitly added constant 
in H is chosen to make H{z, z) zero at these points. 

Now we write down the equations of motion for the instantons 

z = \/2w- 



[l + zz) 



3 

z = -V2w '~' ^ . (5.4) 

We seek a solution which goes from (zi,^j) = {—%■,%) to {zf,Zf) = (z, —i). The two equations 
in ( ^.4|) can be decoupled by exploiting the energy conservation condition H{z, z) = which 
follows from the Hamiltonian nature of the trajectory. This can be written as 

2(^2 + ^2) + 1 + 2^z + zV = 0. (5.5) 

and may be solved to yield z as a function of z and vice versa: 

.^/2z + i .V2z + i 
Z = —I Z = —I (o.d) 

z + V^i z + ^i 
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(Choosing the other solution of the quadratic equation yields instantons running in the 
opposite direction.) Substituting these formulae in the equations of motion yields 



z = —iw{l + z ), z = iw{l + z ). 
These may be integrated by elementary means to yield 



(5.7) 



e^""^ -a 

^ci(t) = -%■ 



+ 



i tanh w(r — tq) 
= —i tanh voir — tL 



0)1 



(5.8) 
(5.9) 



where C = e^"""^", C = e^^'^o. These constants are not independent. Energy conservation 
requires 

c" V2-1 



c + 

It is useful at this point to find the frequencies u, ui and 002. We have 



(5.10) 



2j dzdz 



002 



(5.11) 



2j dz^ 

where the suffix i means that the derivatives are to be evaluated at the initial point. Carrying 
out the algebra, we obtain 



Hence, 



(5.12) 



(5.13) 



We can now evaluate the Wess-Zumino and Solari-Kochetov terms in the tunnelling 
action. We denote these by Jwz and Isk respectively. We begin with Jwz- If we make use 
of Eq. (^), we find 



[ZZ — zz] 



-twiz + z). 



l + zz 

Substituting the explicit forms and performing the integration we get 



(5.14) 
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/wz = -J HC/C) = -2j ln(l + V2). (5.15) 
Now we integrate the Solari-Kochetov term along the instanton trajectory. We find that 

By energy conservation this equals 

^i2 + zz). (5.17) 



Thus, along the instanton, 



w 

')sK-^^i = --^C^-zz) = -iw{z + z). (5.18) 



Hence 



1 r°° r- 

/SK = - / (0SK - ^i)dT = - ln(l + V2). (5.19) 

/ J — CO 

The total tunnelling action is 

/= -(2j + l)ln(l + y2), (5.20) 

The shift 2j ^ 2j + 1 is due to the Solari-Kochetov correction. 

We must now evaluate P. This consists of a product of various factors, all of which are 
to hand. Thus, 

^ ^ (5.21) 



TtN 477 

The factors and (- are found by differentiating the formulas { f).8\j and (|5.9|) and examining 
the limits r ±oo. In this we way we get 

C-C+ = 4^. (5.22) 

Finally, 

2uj 



UJi ~ 00 
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472(3 + 2^2). (5.23) 



Putting these together, we obtain 

P = -il(4 + 3v^)^ = ^v^. (5.24) 

At this point we have almost all that we need to write down the answer for the tunnel 
splitting — except that we need to consider a second instanton. The trajectory (|5.8|) and 



i WM passes close to the north pole of the sphere. By symmetry there must be a second 
instanton which passes near the south pole. This is given by 

Zci = icoth.w{r — tq), Zd = —icothw{T — Tq). (5.25) 

It is obvious by symmetry again that this instanton has exactly the same amplitude as the 
first, so the total amplitude (and thus the splitting) is obtained by simply doubling the 
answer from the first instanton. Hence 

A = lQw( iy^'2i/4e-(2.+i)in(i+v/2)_ (526) 



This agrees with i,|,|l8| [In the last reference put = 1/^2 in Eqs. 4.31-4.34.] In p 
there a numerical comparison which shows that the prefactor is indeed correct. 
For completeness, we note that the average energy is given by E^v = ^{u — ui). 



VI. APPLICATION TO FEs 

The LMG model is of interest primarily because it provides a check of our formalism 
against other well-confirmed calculations. In this section we will calculate the tunnel splitting 
for a family of models that includes a realistic approximation to the molecular magnet Feg. 
The spin-direction dependent energy in Peg is less symmetric than that of the LMG, and 
the relevant Hamiltonian includes an externally imposed magnetic field which serves to pull 
the classical minima off the equator of the unit sphere. It is the experimentally observed 
oscillations in the tunnel splitting as a function of the external field that makes this system 
interesting. The oscillations are a consequence of interference between the two distinct 
instanton trajectories and are accurately reproduced by our calculation. 
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We take as our Hamiltonian 

H = kJl + k2Jl - QI^bHJ,, (6.1) 

with ki > k2 > 0. We define A = k2/ki, He — 2kij/giJLB and 

h = H/H^. (6.2) 

We will express all results in terms of the combinations A and h. It is also convenient to 
define a 1/j corrected field and anisotropy ki by 

h^jh/{j-l), k^^kU-\)l3. (6.3) 

We follow the same steps as in the LMG model. The "classical" Hamiltonian appearing 
in the path integral is 



{z\z) 



(l--.)2-A(.--)2-2/,(l--V) 



{l + zzf 

The energy minima are now at the points 



(6.4) 



z = z = ±^0) (6.5) 

where zq is real and given by 

Zo=[{l-h)/{l + h)Y'^ (6.6) 

In Cartesian coordinates these minima lie in the xz plane — provided we confine ourselves 
to h < 1, which we shall do. In fact, we will assume that 

h < Vl- A. (6.7) 

At the two minima, the energy is 

eo^H{zo,Zo)^-kifh\ (6.8) 
The classical equations of motion are 
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[1 + zz) 



2z{l -zz) + \{z- z){l + z^) + 2hz{l + zz 
-2z{l - zz) + \{z - ^)(1 + z^) + 2hz{l + zz) 



(6.9) 



(1 + zz) 

We wish to solve these subject to the boundary conditions zi = z{—oo) = zq, Zf = z{oo) = 
—zq. Note that Zi = Zi, zj = Zf, so the instanton end points still lie on the real sphere, but 
the rest of the instanton does not. Once again the equations can be decoupled by exploiting 
the fact that energy is conserved along the instanton trajectory. In this case H{z, z) = eg. 
This condition can be written as 



1 - zzy - X{z - zY - 2h{l - z'z') = -h\l + zz 



(6.10) 



and may be solved to give 



y/Xz ±il-h) 



VX±{l + h)z' 

Substituting this in the equation of motion for z, and simplifying, we get 

z = ±y/X{l + h)~ki3{zl-z^). 



(6.11) 



(6.12) 



We will see that to obtain instantons going from Zq to —Zq, we must pick the minus sign in 
this equation. The other sign yields instantons running in the opposite direction. 

It is now elementary to integrate Eq. ( |6.12| ), and use Eq. (|6.11|) to obtain the time 
dependence for both Zci{t) and Zci{t). We find 



^ci(t) = —Zq tanht, 
_ , , tanht + \/l — /i^ 



Here, 



v^A + VI - h^tanht 



t = ut/2, 



(6.13) 
(6.14) 



(6.15) 



and the frequency u is given by 
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(6.16) 



That this is the same uj that follows from Eqs. ( p.l6| ) and (|3.17|) shall be shown shortly. 
It can be seen that our solution corresponds to choosing the minus sign in Eq. ( |6.12| ) as 
asserted above. It is also useful to note that the solution ( |6.13|) and ( |6.14| ) can be rewritten 

as 



Zc\ = —zq tanht, Zd = —zq coth(t + 



(6.17) 



where 



tanh tn 



X 



1/2 



(6.18) 



Equations (|6.12| ) and (|6.11|) possess a second solution 



-Zq coth t, Zc\ = —Zq tanh(t + to 



(6.19) 



Formally, this new trajectory can be obtained from the first by the shift t t + 111/2. 
Alternatively, we could obtain it by switching the expressions for Zd and Zci in Eqs. ( |6.13| ) 
and ( |6.14| ), which corresponds to reflection in the xz plane — a symmetry of the Hamiltonian 
— and then shifting t by — to- 

Again we find the frequencies u, ui and u!2- We note that 



2j 



dzdz 



UJ2 



(6.20) 



2 J dz'^ 

where the suffix i means that the derivatives are to be evaluated at the initial point 1, = z 
Zi. Carrying out the algebra, we obtain 



cji = ^ij(l-/i' + A), 
cu2 = fcij(l-/i'-A). 



(6.21) 
(6.22) 



We now use Eq. ( p.l7| ) to show that uj is given by Eq. ( p.l6| ). The same frequencies are 
found at the final point 1. = z = Zj. 
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We next evaluate and integrate the Wess-Zumino and Solari-Kochetov terms in the tun- 
nelhng action. We denote these by J-wz and I^k- Since the calculations are somewhat 
lengthy, it is best to do the two terms separately. We begin with /wz, considering instanton 
1, i.e., that given by ( |6.13| ) and ( |6.14| ). After some algebra, we obtain 



where 7r2 and tts are polynomials of degree 2 and 3, whose explicit form we do not require. 
What we do need is the differential awz dr. If we make the substitution 

i; = tanht, (6.24) 

and factorize the polynomials tt2 and vra, we obtain 

awzirjdr = - / -dv, (6.25) 

J-l{v-Vi}{V-V2){v-V5) 



where 



.,^-L(i±|)"(_l.^), (0.26) 



^3,4 = ^ , (6.27) 

V, = (6.28) 

1-/^2 



The integral is best done by decomposing the integrand into partial fractions. We find 

- ' +^ ^, (6.29) 



{V - Vi){v - V2){V - V^) V -V5 V ~Vi V - V2 

where 



^=--^. (6.30) 



Thus, 

^ - \ ^, / ^, f 1 - V2 



/wz = j / awzdr = -j In — ) + f3\n[ — - /? In 

J-00 I V— 1— -^5/ \—l—v-[/ 



.-l-v^J \-l-ViJ \-l-V2 
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(6.31) 



The ratio involving is 



1 / . -^1) 

while the 13 terms combine to yield the logarithm of 

1 - ViV2 + {V2 - vi) _ hV\ + Vl^\jl- 

Collecting together the various parts, we have 



(6.32) 



= R2 



(6.33) 



Iwz,i = -jln-Ri + 



jh 



lni?2. 



(6.34) 



We have added another suffix to show that this pertains to instanton 1. 

The next step is to integrate the Solari-Kochetov term. For this we ffist need 0sk- From 
Eqns. 



5SK 



and ( p.2|) we find 



2{i + zzy 



-4 (l - 2 
+4:h{l + zzf 



zz — izz)'^ 



+ 2A ((1 + zzf + 2{z - zf + zz{z - zf) 



(6.35) 



(The reader may verify that as r ^ ±oo, {/)gK — >• uJi- This provides a check on our earlier 
calculation of uoi.) After a little more work, we find, 

ki3 



OSK - ^1 



{l + zzf 



z - z 



i + hy{i + zzY + x{z-zY-4 

2 



(6.36) 



1 + zz 

We have written this expression in such a way that it is convenient to integrate the terms 
on the two lines separately. That is, we write. 



1 

hK = - {4>sK-^^i)dT = Ia + Ib, 

Z J —oo 



(6.37) 



where, 



Ir ■ r , {l + h)\l + zzY + Hz-zy -4: 
2 j-oo (1 + zzY 

Ib = TT^ijA / dr- ^ 



1 + zz 



(6.38) 
(6.39) 
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It is to be understood that the integrands are evaluated along the instanton trajectories. 
We may exploit this fact to simplify the integrand for Ia by using energy conservation and 
so eliminate the term in A. When this is done, we obtain 

= kail + ft) r <i.z(l4±il±l)!£. (6.40) 

J-oo L + ZZ 

The integrals are evahiated in the same way as /wz- With the same change of variables, 
and definitions of Vi to as before, for Ia we get 



VX{l-h) J-1 {v -vi){v -V2) 
(l + h) 



2V1 

(l + h) 



-xJ-i 



dv 



Likewise, for Ib we get 



V — Vi V — V2. 

In ^2- (6.41) 



Is--'-f\ (--^)(-+/) d. (6.42) 



The partial fraction decomposition yields 

{v-l){v + l) 1 /?' /?' 



{v - Vi){v - V2){v - V^) V-V5 V-Vi V-V2 

where 



(6.43) 



(3' = -(1 - A)-^/^ (6.44) 

Hence, 

7B = -^ln^i + ^-^i==ln^2. (6.45) 
Thus the Solari-Kochetov contribution to the action for instanton 1 is 

hK = -\ In ^1 - ^^7^ In R2. (6.46) 

Note that this is 0(l/j) relative to the Wess-Zumino contribution. Adding together the two 
contributions, we obtain the total action 
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/ = -(j + i)ln/?i + 



jh 



InRo 



(6.47) 



In the second term we have used the formula (j — ■^)h = jh. 

We now turn to the prefactor P. In evaluating this, we may ignore differences of order 
i.e., we may replace J by j, h by h, etc. The quantity consists of a product of various 
factors, all of which are already available. Thus, 

(6.48) 



The factors C+ and C- are found by differentiating the formulas ( |6.13| ) and ( |6.14|) and 
examining the limits r ±00. In this we way we get 



C- = -2^0, 



Finally, 



200 ^ ^jKl-h') 
= 4 i 



A(l-/i2) 



ui-uj (1 _ /,2 + ^) _ 2^A(1 - /i2) ^^/r^_^ 
Making use of the identity 



2 • 



2^0 

1 + ^0^ 



we obtain 



P 



4j A^/^(l - /i^)3/2 
' V 1 - /i2 - A ■ 



(6.49) 
(6.50) 



(6.51) 



(6.52) 



(6.53) 



We can now obtain the contribution of instanton 1 to the tunnelling amplitude by sub- 
stituting Eqs. (|6.47| ) and ( p.53| ) in the general formula (|4.10|) . Denoting this quantity by 
Ai, we have 



Ai = 2uJ^f\P\e^-'^/^, 



where the additional factor of e arises from the fact that P < 0. 



(6.54) 
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It remains to obtain the tunnelling amplitude A2 from the second instanton. Because the 
two instantons are related by a complex shift in t, it is apparent that the actions /i^2 (where 
we temporarily add suffixes to distinguish the two) and the prefactors Pi^2 will be given by 
the same analytic expressions. However, the phases to be assigned to the actions and a/P are 
somewhat ambiguous. Unlike the case of a particle moving in one dimension, the prefactor 
in the general formula does not arise as the determinant of a Hermitean quadratic form, 
and there is no unambiguous way for factors of i to get partitioned between the prefactor 
and the exponent. The surest way of fixing the relative phases is to appeal to a physical 
argument. Alternatively, this can be regarded as fixing the signs of the amplitudes and 

For the Feg Hamiltonian ( |6.4|) , let us work in the basis \ j, m) with the standard defini- 
tion of the raising and lowering operators J±, so that the matrix elements (j, m ± 1| J±|j, m) 
are all real. Then the matrix of H is completely real, and since it is Hermitean, all its 
eigenvalues and eigenvectors are also real. Secondly, since Zi = Zq and Zf = —zq are real, the 
states \zij) are real, i.e., all the matrix elements {j,m\zij) are real. Thus the amplitudes 
ttj and af are real. It follows that the amplitude K is real, and so is the one-instanton 
contribution to it, i.e., Ai + A2 is real. Therefore, we must have 

A2 = Al. (6.55) 

Equation (|6.55| ) determines A2, and the energy splitting A completely. However, it is 



still useful to investigate the origin of the phase difference in the actions a little more closely. 
As readers will have noticed already, the integrand in Eq. (|6.25|) is singular at f = f 2 and 



V = V5, since for h < Vl — A, 

vi < -1, -1 < f2 < 1, -1 < f5 < 1. (6.56) 

Correspondingly, both Ri and R2 are negative, and both In Ri and In R2 must be interpreted 
to have an imaginary part of vr modulo an integer multiple of 2tt. The question is what the 
assignment should be for the two instantons. We can see this most easily by examining the 
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difference AJwz = -^wz,2 — -^wz,i- To this end, we note that the WZ one-form may be written 
as a complex one-form in the z plane, 



awz dr 



1 + zziz) 



dz _, , 

'Tz-~'^'^ 



dz = F{z)dz, (6.57) 



with z{z) given by Eq. ( |6.11| ). Thus, Jwz may be written as a ;2-plane contour integral of 



F{z) from zq to —Zq. In fact, apart from a scale factor of 2:05 the substitution ( |6.24D is 



tantamount to changing the integration variable to z, so we see that F{z) has poles at Z0V2 
and Z0V5 (the one at zqVi does not matter). The two instantons go around these poles in 
opposite senses, so A/wz is given by integrating F{z) along a closed contour from zq to —zq 
and back to zq: 

AIyvz = f F{z)dz. (6.58) 

The residues at the poles can be read off the partial fraction decomposition ( |6.29|) , yielding 

2jh7r 

%z,2 - Iwz,i = 2j7r - (6-59) 

This is precisely what we would obtain from Eq. ( |6.55| ), for that would have us assign iivr 
for lni?i (and lni?2) for the two instantons. 
The energy splitting is given by 

A = Ai + a;. (6.60) 

To compare with previous results, it is useful to rewrite this as follows. Consider the real 
part of the action. 

To = -Re / = (j + i) In - -^^= In |i?2|. (6.61) 

The ratios Ri and R2 are defined in terms of the field h. If we write h = h + 0(1/ j), and 
expand in powers of we discover that 

To = (j + ^) In \R^\ - -^i= In \R2\ + Oir'), (6.62) 
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where Ri is obtained from Ri by simply deleting the tildes above the /I's. Note that the 
corrections are of 0(l/j), not 0(1). These are beyond the accuracy to which we are working, 
so we simply drop them henceforth. 

Thus, the complete expression for the splitting is 



A = i/-cuF^/V^° cos A. 

TT 



(6.63) 



We give the expressions for F, Fq and A for ready reference: 

.AV2(1_/,2)3/2 



1 - /l2 - A 



To = {3 + -)\n 

A = Im / 

2 

= JTT 1 



jh 



In 



1- A)(l-/i2) -/i^A 



(6.64) 
(6.65) 



h 



(6.66) 



Our answer for A is identical to that found by means of the discrete WKB method in 
||T9| [see Eqs. (5.1-5.5)]. Naturally, the points at the which the tunnel splitting vanishes are 
the same too. 

The nontrivial aspect of this calculation is that there are 1/ j corrections in the quenching 
condition. If we simply take the energy expectation H{z,z) = {z\H\z) / {z\z) in the Wess- 
Zumino term, we have the problem that the anisotropy and field terms scale with j differently 
if 1/j corrections are included. This is how the quenching condition was found in [Q, but the 
1/j corrections were never considered, so it was some what serendipitous that the condition 
that was stated turned out to be rigorously correct. By including the SK correction, this 
deficiency is now repaired. 



VII. CONCLUSION 

We have seen that, once the extra-phase contribution is included, the coherent-state path 
integral for spin provides an accurate and effective tool for calculating tunnel splitting. In 
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particular the Weyl shift j ^ j + | appears automatically. It must therefore be possible to 
put the spin coherent-state path integral on the same sound mathematical footing as the 
conventional Feynman integral. 
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IX. APPENDIX 



Here we derive equation (|4.6|). We first apply an SU{2) rotation to 



^^initiaK^;, Z) = ^\ UJi{z - Zi){z - Z*) + \uj'2,{z - Zif + \ujI{z - Z*f (9.1) 

\L -\- Z^Zi) 

in order to place Zi^ "Zi at the origin, and to make the coefficient uj2 real. The result is 

H{z, z) = 2j [uizz + \uj2Z^ + \uj2z'^\ . (9.2) 

In the semiclassical limit, 2j ^ 1, we may ignore the curvature of the phase space and, 
after rescaling ^/2J z z to account for the difference in the coefficient in the kinetic terms, 
identify H{z, z) with the coherent state classical Hamiltonian for the squeezed harmonic 
oscillator 

H = uia)a + \u2{a^'^ + 0^). (9.3) 
The Bogoliubov transformation 

h = cosh 9 a + sinh 6 a) 

fel' = sinh 6' a + cosh 6^ (9.4) 

reduces the Hamiltonian 
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^squeezed = ^ cosh 20 (ci^ a + |) + sinh 2e{a)'^ + a^) (9.5) 

to 

Squeezed = ^ {b^b + i) , (9.6) 

and so we identify 

Q = uj = \Juj'l — UJ2, 
Q cosh 26 = uji, 

fisinh2^ = cj2. (9.7) 
The eigenvalues of H are therefore 

E„ = o^(n+ i) - icui. (9.8) 
The operators a^a, and a^^ generate the Lie algebra su(l,l). Therefore either the flat 



phase-space coherent state path integral or standard sm(1, 1) disentangling methods |2Q| . |21 
can be used to derive 

(C/|e-^^|0) = D--2 exp {D-' {Cfd - \ sinh 20 sinh o.T(c' + C- )) } e"^'^^^, (9.9) 

where 

D = e"^^ cosh^ e - e-'^^ sinh^ 9, (9. 10) 
and the harmonic oscillator coherent states \() are defined by 

10= exp Cat |0), a|0) = 0. (9.11) 
In the large— T limit, and with Q and C/ both at the origin, this reduces to 



(0|e-^^|0) ^ (cosh0)-ie-^(--"^)^ = J ^^e-'^^^^^^^ . (9.12) 

^ UJ -\- UJi 

We now rotate back to the original Zi. Taking note of the transformation properties of the 

|2;)'s, we get 



= (1 + z.^.)'^/^^e-i(--)^, (9.13) 
y UJ -\- uji 



as claimed. 
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